Influence of gas pressure on high-order harmonic generation of Ar and Ne 
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We study the effect of gas pressure on the generation of high-order harmonics where harmonics due 
to individual atoms are calculated using the recently developed quantitative rescattering theory, and 
the propagation of the laser and harmonics in the medium is calculated by solving the Maxwell's 
wave equation. We illustrate that the simulated spectra are very sensitive to the laser focusing 
conditions at high laser intensity and high pressure since the fundamental laser field is severely 
reshaped during the propagation. By comparing the simulated results with several experiments 
we show that the pressure dependence can be qualitatively explained. The lack of quantitative 
agreement is tentatively attributed to the failure of the complete knowledge of the experimental 
conditions. 
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I. INTRODUCTION 



In the last two decades, high-order harmonic genera- 
tion (HHG) has been widely studied for its potential as 
a short-wavelength tabletop light source [iry], or as an 
ultrashort attosecond pulse or pulse train [J,Q. The har- 
monic emission in a gas medium is well understood quali- 
tatively. Harmonics are emitted when atoms or molecules 
are exposed to an intense infrared laser field. These lights 
travel collinearly with the input laser and interact with 
the generating medium. Thus a full description of HHG 
involves a quantum treatment of the harmonic emission 
from a single atom or molecule, together with the propa- 
gation of the harmonics and the nonlinear interaction of 
the laser light with the medium. Recently, Jin et al. [Q- 
1^ have combined the well-established propagation theory 
with the quantitative rescattering (QRS) theory |10l - KL3 | 
for individual atoms and molecules, to provide a quan- 
titative description of HHG generated from atomic and 
molecular targets by an intense laser pulse. Its success 
has been shown by comparing with the recent experi- 
mental measurements _£., J_, 2\ ■ These applications have 
focused on the region where the fundamental laser field 
is not severely modified in the medium. Such studies are 
useful for using HHG to probe the structure of molecules. 
For applications, it is desirable to achieve highest num- 
ber of photons by increasing the laser intensity and the 
gas pressure. If the HHG from individual atoms are fully 
phase-matched, then the number of photons will increase 
quadratically with the medium pressure. Unfortunately, 
the nonlinear interaction of lasers with the gas medium is 
complicated. The intense laser light and the harmonics 
can be dispersed and absorbed by the medium. The in- 
tense laser light also ionizes the atoms (or molecules) and 
generates free electrons, thus changing the optical prop- 
erties of the gas medium. These nonlinear interactions 
of the laser light and the harmonics with the medium 
can be calculated by solving the proper Maxwell's wave 



equations, with the induced dipole from each atom or 
molecule by the laser as the source of the harmonics. A 
detailed theory of harmonic generation with the inclu- 
sion of macroscopic propagation where single-atom HHG 
emissions are calculated with the QRS theory has been 
given in Ref. Q. 

In this paper our goal is to extend this theory to the 
harmonic generation for laser intensities near and above 
the critical intensity for a given target. In particular, we 
want to investigate the effect of pressure on the yield of 
HHG, as the intensity is varied or when the focusing con- 
ditions are changed. At high intensities, HHG is a highly 
nonlinear process, and the harmonic yields are very sen- 
sitive to all the parameters of the experiment. These 
parameters, unfortunately, are not always well specified 
in an experiment. We perform simulations by varying 
some parameters in order to achieve optimal agreement 
with the reported experimental data. We also check the 
effects of various parameters of the medium to illustrate 
how these individual parameters change the simulated 
HHG spectra. These kinds of studies have been carried 
out previously, sometimes connecting with experimental 
observations. However, in these calculations HHG from 
each atom was often calculated using the strong field ap- 
proximation (SFA), or the so-called Lewenstein model 
(l3| . It is well known that the SFA does not predict 
correct HHG spectra from individual atoms. By using 
QRS, the atomic dipoles induced by the laser are accu- 
rately calculated. These induced dipoles are then fed 
into the Maxwell's propagation equations. We perform 
simulations to check to what extent some of the earlier 
experimental HHG spectra can be simulated with the 
present theory. 

In Section [Hi we briefly summarize how the calcula- 
tions are done. We then simulate the HHG spectra of 
Ar for 6-fs pulses at quite high intensity reported in Rcf. 
|14| . We also examine the effect of gas pressure on the be- 
havior of Cooper minimum at high intensities. For higher 



pressure the laser field is severely modified as it propa- 
gates through the medium. We show that the Cooper 
minimum eventually disappears at pressure above about 
50 Torr. The disappearance of the Cooper minimum is 
entirely due to the medium propagation effect [15]. Next 
we examine the pressure effect on the harmonics, and it 
is established in general that there is an optimum pres- 
sure for each harmonic beyond which the harmonic yield 
drops with increasing pressure. Since the propagation 
of harmonics is controlled by the optical properties of 
the medium, we check if there is a single factor that 
would mostly affect the harmonic generation. We also 
take a close look at the dispersion and absorption coef- 
ficients used in the simulation. These tabulated data, 
which were obtained from theoretical calculations, have 
not been carefully calibrated and should be used with 
caution. Finally we examine the propagation of HHG 
generated in Ne at high intensity and high gas pressure 
to compare with the result reported in a recent experi- 
ment [lal ■ We have found that in this case HHG spectra 
are very sensitive to the experimental parameters such 
that simulation can only achieve qualitative agreement. 
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FIG. 1: Comparison of theoretical and experimental HHG 
spectra of Ar generated by a 760-nm and 6-fs laser. Ex- 
perimental data are from Ref. [lj|. For simulation- 1 and 
simulation-2 the laser peak intensity at focus is 4.5 and 
4.0x10^* W/cm^, the length of gas jet is 1 and 2 mm, re- 
spectively. The other laser parameters are given in the text. 



II. THEORETICAL METHOD 



III. RESULTS AND DISCUSSIONS 



The detailed description of the theoretical method used 
in the calculation can be found in Ref. Q. Briefly, we 
consider both the propagation of the fundamental and 
the harmonic fields in an ionizing medium. For the fun- 
damental field, we include dispersion, absorption, Kerr, 
and plasma effects. For the harmonics, the dispersion 
and absorption due to neutral atoms are included. The 
induced dipoles for single atoms are obtained using the 
QRS theory [ll|: 



i?QI^S(^)^p^SFA(^)^QRS(^)^ 



(1) 



Here PF^^^(cj) is the wave packet calculated from the 
Lewenstein model [l3|, normalized by the tunneling 
ionization rate calculated using the Ammosov-Delone- 
Krainov (ADK) theory [17|, l18|. The transition dipole 
moment d'^^^{uj) is obtained by treating the real atom 
in a single-active electron approximation. For the argon 
target used in the present work, d'^^^{uj) is calculated by 
using the model potential given by Miiller [19|. For Ne, 
the model potential takes the form given in Refs. |8l.[l8|. 
The resulting induced dipoles for hundreds of different 
peak intensities are then fed into Maxwell's wave equa- 
tions. We assume that the laser beam at the entrance 
of the gas jet has the Gaussian shape both in time and 
space. The harmonics emitted at the exit face of the gas 
jet may propagate further in free space until they are 
detected. Thus the harmonic spectra observed experi- 
mentally also depend on the detecting conditions 01 . In 
this paper, except for one case to be described below, all 
the harmonic signals (i.e., the total signal) from the exit 
of the medium are collected. 



A. Macroscopic HHG spectra of Ar: Theory vs 
experiment 

We first simulate experimental HHG spectra of Ar gen- 
erated by a 760-nm and 6-fs (FWHM) laser pulse re- 
ported in [Tjj . Experimentally, the gas jet was placed 2 
mm after the laser focus. The confocal parameter of the 
laser beam was 10 mm (the beam waist was 34.8 /xm if a 
gaussian beam is assumed). The laser intensity at the fo- 
cus was estimated to be 6x 10^"* W/cm^. In simulation-1, 
we use the gas-jet width of 1 mm, and the gas pressure is 
40 Torr. We find that a laser peak intensity of 4.5x10"^^ 
W/cm^ would give better agreement with the spectra in 
the cutoff region. The calculations have been carried out 
for different carrier-envelope phases (CEPs) since a short 
6-fs laser pulse was applied in the experiment. The theo- 
retical spectra are finally averaged over all the CEPs. To 
compare with experiment, the HHG spectra are normal- 
ized at harmonic 27 (H27). We can see that the main fea- 
tures in the experimental spectra are well reproduced by 
simulation-1. However, the overall slope of the harmonic 
yields is not the same. The simulated harmonic yields are 
too low for lower orders, but too high for higher orders. 
But this result is already a great improvement over the 
theoretical simulation reported in Fig. 4(c) of Ref. |14| . 
Since HHG spectra are sensitive to experimental condi- 
tions, we thus make another simulation where the gas-jet 
length is extended to 2 mm and the laser intensity is 
reduced to 4.0x10^'* W/cm^. In this simulation-2, other 
parameters are the same as those in simulation-1. We can 
see that simulation-2 gives better agreement for the low 
harmonics up to H27, but the slope for higher harmonics 



is still not reproduced correctly. We have made further 
adjustment of the laser parameters within the capability 
of our present code but were unable to improve agree- 
ment. At present our code is limited to a Gaussian (spa- 
tial) input pulse. In the experiment of Altucci et al. [14| . 
the laser was coupled into an argon-filled capillary con- 
sisting of a 60-cm-long fiber whose spatial profile would 
be the Bessel function. Thus the emerged 6-fs pulse may 
be closer to a truncated Bessel beam. It remains to be 
seen if the discrepancy can be removed if a truncated 
Bessel beam is assumed for the input pulse |20l - [23 |. 

Below we aim at investigating the gas pressure effects 
on the HHG spectra. 



B. Influence of the gas pressure on Cooper 
minimum and spectral modulation 

In Figs.UJa) andj^fb) we show the macroscopic HHG 
spectra of Ar for four different gas pressures: 10, 30, 50, 
and 80 Torr, at a laser intensity of 5x10^^ W/cm^ at 
focus for GEP=0 and GEP averaged, respectively. The 
other parameters are the same as those in simulation- 
1 of Fig. [1] Since a few-cycle pulse is used, the har- 
monic spectra show a strong GEP dependence in the cut- 
off region. After GEP average the spectra become very 
smooth. There are two additional features in Figs. [DJa) 
and^b) that warrant further discussions. 

The first is about the well-known Gooper minimum 
(GM) of the HHG in Ar. Ahhough GM in harmonic 
spectra is a signature of the electronic structure of the Ar 
atom [23, [2J] , previous experimental studies have shown 
that the "visibility" of GM is sensitive to experimental 
conditions. For example, the position of the gas jet 25 1 
and the gas pressure [l^ could wash out the GM seen 
in the harmonic spectra. It has also been pointed out 
by Jin et al. [3] that the GM also depends on how the 
HHG spectra are measured, i.e., in the near or the far 
fields, with or without the slit. The present simulations 
show a similar GM dependence on gas pressure, which has 
been observed in Ref . [l5| , although the experiment was 
performed under different conditions (laser duration was 
^^50 fs, and the laser intensity was ~2.8xl0^^ W/cm^). 
In Figs. [UJa) and[2][b), there appears an obvious min- 
imum around 31st harmonic (~51 eV) below 30 Torr. 
When the pressure is increased to 50 Torr, the mini- 
mum becomes barely visible, and then it is totally washed 
out when the pressure is increased to 80 Torr. Accord- 
ing to the QRS [y, 0] , these variations are attributed to 
the change of the macroscopic wave packet (MWP). The 
QRS states thaM;he macroscopic HHG spectrum can be 
expressed as 



Si,{uj)ocuj^\W{oj)\^\d{uj)\' 



(2) 



where W{ijj) is the MWP, which reflects the effect of laser 
and experimental conditions, and d{uj) is the photore- 
combination (PR) transition dipole moment, which is a 
property of the target only. In Fig. ^c) we show the 



detailed comparison of the MWPs at four gas pressures 
for GEP=0. It shows that the MWP (|Vl^(w)|) varies 
with the gas pressure. Around photon energy of 50 eV 
(--H31), the MWP is very flat at 10 Torr. With increas- 
ing gas pressure it becomes more steeper. Meanwhile, 
we also show the PR transition dipole moment |d(w)| for 
reference. It is concluded that the rapid change of the 
MWP around 50 eV causes the disappearance of the GM 
in the HHG spectra. 

The second feature in Fig. [2] is the spectral modula- 
tion. In Fig. [UJa), for harmonics above the GM, the 
spectra modulation is large at 10 Torr. The spectra are 
less modulated at 30 and 50 Torr, and become almost 
continuous when the gas pressure is increased to 80 Torr. 
These continuous HHG spectra could be used to produce 
isolated attosecond pulses (lAPs) with a GEP-stabilized 
laser. Note that the continuum structure in HHG spectra 
of Xe has been discussed by Jin et al. [2a] recently [see 
Fig. 1 in that paper]. They found that reshaping (blue 
shift and defocusing) of the fundamental laser field was 
responsible for these phenomena, and they also showed 
a method to produce an lAP by spatial filtering in the 
far field. To check this mechanism in our situation, we 
show the evolution of the on-axis electric field at the en- 
trance and exit of the gas jet for 10 Torr and 80 Torr in 
Fig. md). Note that the time is defined in the moving 
coordinate frame [7|, |26j. The ionization probability is 
very high (about 61% at the end of a laser pulse with 
peak intensity of 4.63x10^^ W/cm^). The electric field 
has a good Gaussian form at the entrance. In the lead- 
ing edge of laser pulse, where the ionization probability 
is very small as seen in the inset of Fig. [2jd), the electric 
field at the exit face has a small shift with respect to the 
one at the entrance. This shift is originally from the ge- 
ometric phase due to the tightly focused laser beam. In 
the falling edge, there is an obvious blue shift, and the 
electric field at 80 Torr is much reduced than that at 10 
Torr. We can conclude that the fundamental laser field 
is reshaped. Gombined with effects due to GEP average, 
our simulations in Fig. [UJb) show that spectral modula- 
tion above the GM increases from 10 Torr to 30 Torr, 
and then decreases with increasing gas pressure until the 
HHG spectra become continuous at 80 Torr. Similar phe- 
nomena have been observed in Zheng et al. [271 ] for a 
phase-stabilized few-cycle laser pulse. 



C. Gas pressure dependence of HHG conversion 
efficiency at fixed laser intensity and gas-jet length 

1. The optimal HHG yield vs gas pressure 

The macroscopic harmonic generation is strongly in- 
fiuenced by many factors collectively. In other words, 
the conversion efficiency cannot be improved by just in- 
creasing the laser intensity, atomic density, and medium 
length independently. There are three physical effects 
that limit harmonic generation: absorption, defocusing. 
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FIG. 2: Macroscopic HHG spectra of Ar at the gas pressures of 10, 30, 50, and 80 Torr for (a) CEP=0 and (b) CEP averaged. 
(c) The normahzed macroscopic wave packet (MWP) at four different pressures for CEP=0 and transition diploe moment 
|d(aj)|. (d) Evolution of the on-axis electric field: at entrance and exit. The inset shows the ionization probability vs time at 
entrance. The laser peak intensity at focus is 5xlO^*W/cm'^. The other parameters are the same as those in simulation-1 in 
Fig. [II 



and dephasing [28|. Many experimental and theoretical 
studies that investigated harmonics efficiency (or pho- 
ton flux) have been carried out, such as absorption by 
Schniirer et al. [23, |30| , defocusing by Altucci et al. [3l| 
and Dachraoui et al. [la], and dephasing in a hollow 
waveguide in Refs. [33.I33J. 

Figure |3Ka) shows the dependence of harmonic inten- 
sity on the gas pressure for different plateau harmonic 
orders between H13 and H21 for CEP=0 at the inten- 
sity of 2x10^'* W/cm^. The gas jet is 1-mm long, with 
other parameters the same as those in simulation-1 of 
Fig. [T] With the increase of gas pressure, the yield for 
each harmonic order increases first to reach a maximum 
and finally drops with even higher pressure. The maxi- 
mum yield shifts to a higher pressure for higher harmonic 
order, especially in the region of low harmonic orders. 
These have been demonstrated experimentally in an Ar 
gas-filled hollow waveguide [30| and in a gas cell [3J| . We 
show a similar figure for CEP=7r/2 in Fig. [3Kb). 

The above results are most easily understood from a 
one-dimensional model given by Constant et al. [35|. 



For the qth. harmonic order, the number of photons is 
proportional to [Eq. 1 of Ref. [35| ] 
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where Aq{z) (units: C-m) is the amplitude of the atomic 
response. Here Lcoh=7i'/Afcg is the coherence length, 
where Afc^ is the wave vector mismatch between the fun- 
damental and the generated harmonic field, Labs=l/cp is 
the absorption length, where a and p are the photoion- 
ization cross section and the density of the generating 
gas, respectively. According to this model, if we assume 
Aq{z) to be independent of z (this is true for a loosely fo- 
cused laser beam) , then the harmonic yield is determined 
by the three parameters: the medium length -Lmed, co- 
herence length Lcohj and absorption length Labs- The 
phase mismatch Afcq is determined by different disper- 
sion terms: atomic, electronic, geometric dispersion (the 




30 40 

pressure (Torr) 




20 30 40 50 

pressure (Torr) 



600 



500 



S 400 



300 



200 



100 




10 20 30 40 50 60 

pressure (Torr) 




10 12 14 16 18 20 22 24 26 28 30 32 34 
Harmonic order 



FIG. 3: The harmonic intensity of Ar as a function of gas pressure for different orders from 13rd to 21st in the plateau for 
(a) CEP=0 and (b) CEP=0.57r obtained by QRS. The laser peak intensity at focus is 2x10^* W/cm^, the other parameters 
are the same as those in simulation-1 in Fig. [T] (c) The same as (a) but obtained by SFA. (d) Comparison of HHG spectra 
obtained by QRS and SFA for a pressure of 30 Torr and CEP=0. 



so-called Gouy phase for a Gaussian beam), and by the 
gradient of the atomic dipole phase [la|. In Fig. [3l the 
medium length is fixed. An increase of pressure makes 
the absorption length smaller. According to Fig. 1 of ^35j , 
the harmonic yield would reach saturation, the larger the 
phase mismatch, the smaller the number of photons gen- 
erated for the harmonic. The results from Fig. [3] which 
are simulated from the 3D model, show that as the pres- 
sure increases beyond the optimal pressure for the gener- 
ation of a given harmonic, continuing increase of pressure 
makes the harmonic yield drop. In fact, since the coher- 
ent length varies spatially for each harmonic, the spa- 
tial distribution of the harmonic signal also changes with 
pressure. In fact many factors will influence coherence 
length or absorption length through their dependence on 
gas pressure. Furthermore, the coherence length and ab- 
sorption length also depend on harmonic order. All of 
these can affect the change of optimal pressure with har- 
monic order. If the single-atom harmonics are calculated 
using SFA, the pressure dependence of each single har- 
monic will be qualitatively correct, see Fig. ^c). On the 
other hand, the relative intensities between the harmon- 
ics would not be correct since SFA does not use correct 
photoabsorption cross sections. Thus the HHG spectra 



predicted using SFA and QRS are different, see Fig.IHJd). 



2. Dependence of harmonic spectra on the optical 
properties of the gas medtum 

In the calculations above, few-cycle laser pulses are 
applied and laser intensity is low such that free electron 
density in the gas medium is very small. We have checked 
that effects due to free electrons (defocusing and phase 
mismatch) in these calculations are negligible. In the sim- 
ulation we include the dispersion ((5i), Kerr nonlinearity 
(772) of the fundamental field, the dispersion (i5h) and ab- 
sorption (/3h) of the harmonic fields in the simulation [3]. 
All of these factors are pressure dependent. Note that 
5i, 772, and (Sh all contribute to phase mismatch A/c^; in 
other words, they all affect the coherence length Lcoh- In 
the meanwhile, /3h leads to pressure-dependent Labs- To 
see how each term contributes to the HHG yield, we re- 
move each term successively from our model. The results 
(envelope of HHG spectra only) are shown in Figs. IH^a) 
andHJb) for two pressures, 10 Torr and 50 Torr, respec- 
tively. At 10 Torr, i5i, 772, and Jh are small, their effects 
can be neglected, but the absorption term /3h is impor- 
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FIG. 4: The influence of dispersion of the fundamental (Si) 
and harmonic (5h) field, absorption of harmonics (/3h), and 
Kerr effect at (a) 10 Torr and (b) 50 Torr. The parameters 
are the same as Fig. |3Ja). (c) The atomic scattering factors 
/i and /2 are from NIST [363 for Ar. 



tant, especially for low harmonics. [Note that /2 is large 
below 25 eV, see Fig. Hlje)]. When the gas pressure is 
increased to 50 Torr, the absolute values of 5i, 772, and 
(5h all increase by five times. Their effects become signif- 
icant. In other words, with increasing gas pressure, the 
change of coherence length cannot be attributed to one 
single factor only. 

For the absorption effect, we first calculate the absorp- 
tion length Labs for II15. They are about 2.4, 0.8, and 
0.5 mm at the gas pressures of 10, 30, and 50 Torr, re- 
spectively. According to Eq. ^ , only when Lmod > ^abs 
the absorption effect would become important. This has 
been verified in our calculation by comparing HI 5 at two 
gas pressures without /3h in the model. The absorption 
effect becomes significant for all harmonics in the plateau 
at high pressure (50 Torr). It is important only for low 
harmonics at low pressure (10 Torr). 

We also show atomic scattering factors /i and /2 in 
Fig. SKc), they are related to ^h and /3h respectively Q. 




30 40 50 60 
Photon energy (eV) 



80 




20 30 40 50 60 70 80 
Photon energy (eV) 
20 25 30 35 40 45 50 55 




15 20 25 30 

Harmonic order 



35 



FIG. 5: (a) Comparison of atomic scattering factors /i from 
Henke [33 and NIST ,36] for Ar. (b) Comparison of photoion- 
ization cross section between experiment [3y| and theory from 
Henke [3^ and NIST [s^ by using the formula a = 2roA/2. 
(c) Comparison of harmonic spectra calculated by using dif- 
ferent /i, /2 from [331 and [3y|, the laser peak intensity is 
2.0x10^* W/cm^, and gas pressure is 50 Torr. 



We can see that the change of harmonics with or without 
Sh (or /3h) follows the energy dependence of /i (or /2). 
From Fig. 21 we conclude that there is not a single dom- 
inant factor that determines the pressure dependence on 
the HHG observed in Figs. |3Ka) and[31[b). Since these 
factors are wavelength dependent, the effects vary with 
the harmonic order. 



3. Dispersion and absorption data from different sources 

In Fig. m we have used the atomic scattering factors 
/i and /2 from NIST [3g| to calculate the dispersion (i5h) 
and absorption (/3h) effects in the propagation equation. 
An alternative set of atomic scattering factors are given 
by Henke ]37\. These two sets of data are compared in 
Figs. [5ja) andjSJb). Since /2 is related to atomic pho- 
toionization cross section through cr ^ 2roA/2, where rg 
is the classical electron radius, we can obtain theoretical 
a from these two sources to compare with experimen- 
tal photoionization cross sections tabulated in [38|. We 
note that the NIST data underestimate the absorption 
cross sections quite significantly. This would have the 
effect of overestimating the simulated HHG yields. Us- 
ing these two sets of atomic scattering factors to simu- 
late the HHG spectra, as shown in Fig. [5{c), there are 
significant differences between the two calculations for 
the lower harmonics. Such differences have not been ad- 
dressed in the literature as far as we know. However, the 
data set from Henke [37] did not extend /i below 29 eV. 
In fact, it drops precipitously below about 40 eV. These 
atomic scattering factors were calculated originally for 
x-ray energy region. They are perhaps not very reliable 
in the photon energies considered here. The accuracy of 
scattering factors in this energy region probably should 
be used with caution. 



D. Pressure dependence of HHG at different laser 
intensities and gas-jet lengths 

Figures[6ja) and[6l^c) compare the pressure dependence 
of H15 and H21 at three different intensities in a gas jet 
with the length of 1 mm. For each intensity, the pressure 
dependence is similar to what we have seen in Fig. [31 As 
the intensity is increased, we note the highest harmonic 
yield occurs at a higher pressure. Comparing the two 
harmonics, the optimum harmonic yield also moves to a 
higher pressure for the higher harmonic. For intensity of 
4.0x10^^ W/cm^, H21 achieves optimal yield at a pres- 
sure close to 100 Torr. Because higher intensity results 
in higher free electron density, consequently, a higher gas 
pressure is needed to compensate for the phase mismatch 
due to the free electrons. This means that if we want to 
obtain high intensity for the high harmonics, we should 
use a high laser intensity and a high gas pressure. This 
result is consistent with the generation of water window 
harmonics using mid-infrared lasers reported recently by 
Popmintchev et al. J-, _2] , where gas pressure as high as a 
few atmospheres were used, together with a higher laser 
intensity. 

In Figs.ini^b) andini^d) we show the harmonic intensities 
for HI 5 and H21 at several different lengths of the gas 
jet (0.5, 1, 2, and 3 mm) using laser intensity of 2x10^^ 
W/cm^ at focus (for each length, the center of gas jet is 
always placed at 2 mm after laser focus). Below 20 Torr, 
the harmonic intensity increases rapidly by increasing the 



interaction length and gas pressure. When we further in- 
crease the pressure to a higher value, on the contrary, the 
harmonic intensity drops with the medium length. The 
optimal gas pressure increases with decreasing interac- 
tion length, because the absorption length is decreased 
with increasing gas pressure. The optimized harmonic in- 
tensity with a long gas jet is higher than that of a short 
one. Thus we can obtain a high harmonic intensity using 
a long gas medium with low gas pressure. These results 
are very similar to those observed by Tamaki et al. [39[ 
using a 30-fs laser pulse. 



E. Macroscopic HHG spectra of Ne: Theory vs 
experiment 

To test how the QRS model works for other atomic 
targets, we perform a simulation of HHG spectra for Ne 
reported by Dachraoui et al. [1^1 recently. Figure [Tlja) 
shows the comparison of HHG s pec tra between theory 
and experiment. The experiment [16| was performed us- 
ing a laser with pulse duration of 45 fs, central wave- 
length of 805 nm. The laser focusing position was in the 
middle of the 3-mm long gas cell and the spot size was 
estimated to be ~110 /im (the beam waist at focus was 
55 /im). The peak vacuum intensity was 1.2(±0.4)xl0^^ 
W/cm^. Three different pressures: 38, 83, and 158 Torr 
were used. In the simulation the beam waist is 65 /im, 
the laser peak intensity is 8x10^^ W/cm^ (near the crit- 
ical intensity of Ne, which is ~8.5xl0^^ W/cm^), the 
other parameters are the same as those given in the ex- 
periment. In Fig. [Tl^a) simulated spectra at the three 
pressures are shown, together with the envelope of the 
calculated spectra. The general trend of the spectra vs 
gas pressure agrees reasonably with the experiment, see 
inset. Under present conditions the ionization probabil- 
ity is very high (^75% at the end of laser pulse with peak 
intensity of 8x10^'' W/cm^ according to ADK tunneling 
model [i3])- The fundamental field is strongly reshaped 
during the propagation in the medium and has a strong 
dependence on pressure. By increasing the gas pressure, 
the harmonic cutoff position shifts to lower orders. 

There are still large discrepancies between the experi- 
ment and the simulation, especially the cutoff positions. 
The cutoff positions in the theory are about 10 har- 
monic orders higher than those in the experiment. In 
our simulations, we find that many experimental param- 
eters would affect the shape and the cutoff of the spectra 
significantly. In Fig. [TJ^b) we show the infiuence of beam 
waist and laser intensity on the spectra at a pressure of 83 
Torr, the other parameters are fixed as those in Fig.[71[a). 
We can see that the shape and cutoff position obtained 
under different conditions are very different. The change 
of cutoff position for the beam waist from 65 to 45 /xm is 
about 10 harmonic orders. A slight change of the laser in- 
tensity from 8 to 7xl0^''W/cm'^ would modify the shape 
of the HHG spectra drastically. Since many factors men- 
tioned above in the experiment are not well determined. 
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FIG. 6: The intensity of harmonics 15th and 21st as a function of gas pressure for (a) and (c): different laser intensities, and 
(b) and (d): different gas jet lengths. In (a) and (c) the gas jet is l-mm long. In (b) and (d) the laser intensity at focus is 
2.0x10^* W/cm^. The other parameters are the same as simulation-1 in Fig. [T] 



we have not been able to simulate the observed spectra. 
We also comment that the experimental HHG spectra 
depend critically on how they are measured [7|. For ex- 
ample, if a slit of width of 500 [iva. is placed at 50 cm 
behind the gas jet, then only harmonics emitted close to 
the axis will be measured. The second inset in Fig. El^a) 
shows how the envelopes of the harmonic spectra at the 
three pressures are modified. For higher pressure, there 
are more free electrons generated and thus the fundamen- 
tal laser is more defocused as the pressure is increased. 
The inset shows that the relative yield at the higher pres- 
sure is reduced, indicating that HHG generated at higher 
pressure is more divergent. 



IV. SUMMARY AND OUTLOOK 

In summary, we have used the recently developed 
quantitative rescattering (QRS) theory for the generation 
of high-order harmonics of individual atoms with the in- 
clusion of the propagation of the laser and the harmonics 
in the gas medium to obtain macroscopic HHG spectra 
that can be compared directly with experimental obser- 
vations. We have studied the effect of the gas pressure on 
the yields of the harmonics. The gas pressure affects har- 



monics through frequency-dependent properties of the 
medium such as dispersion, absorption, nonlinear Kerr 
effect, and plasma density, as well as the phase-matching 
conditions. In the meanwhile, the fundamental laser field 
is modified as it propagates through the medium. The 
modification is especially severe when the incident laser 
reaches the intensity near or above the critical intensity 
for a given target. From the simulations, we have found 
that HHG spectra are very sensitive to the laser focus- 
ing conditions. By starting with the QRS theory instead 
of the SFA for the single-atom response, we have been 
able to obtain HHG spectra that are much closer to ex- 
perimental observations, but we are still unable to obtain 
HHG spectra that are in perfect agreement with the data. 
This is in contrast with our recent positive results [a, S Q 
where HHG of atoms and molecules were generated at 
lower intensities. We have attributed the lack of good 
agreement to the strong dependence of the HHG spectra 
on the laser parameters when experiments are carried out 
at high intensities and that these parameters are difficult 
to specify accurately in a given experiment. In the fu- 
ture, it is desirable that experiments report HHG spectra 
taken at various focusing conditions and gas pressures in 
order to establish if our current theory of HHG is under 
solid ground in the high-intensity and high-gas-pressure 
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FIG. 7: (a) The simulation of total HHG spectra of Ne for 
three gas pressures, 38, 83, and 158 Torr, generated by an 
805-nm and 45-fs laser field. The laser peak intensity at fo- 
cus is 8x10^* W/cm'^. The first inset shows the envelope of 
experimental spectra [lg|. The second inset shows the theo- 
retical spectra if the harmonics are taken with a slit of width 
of 500 pm placed 50 cm after the gas jet. (b) The influence of 
beam waist and laser intensity on the simulated spectra for a 
pressure of 83 Torr. Other parameters are given in the text. 



regime. This regime is important in order to generate 
high-energy harmonics near and above the water window 
using mid-infrared lasers, as well as the possible genera- 
tion of few-cycle single attosecond pulses in the soft X-ray 
region [lM|. 
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